CalrOs post-perovskite, a j — 1/2 quasi-one-dimensional antiferromagnet 
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The 5d 5 iridate Calr03 is isostructural with the post-perovskite phase of MgSiC>3, recently shown 
to occur under extreme pressure in the lower Earth's mantle. It therefore serves as an analogue of 
post-perovskite MgSiOa for a wide variety of measurements at ambient conditions or achievable with 
conventional multianvile pressure modules. By multireference configuration-interaction calculations 
we here provide essential information on the chemical bonding and magnetic interactions in CalrOa. 
We predict a large antiferromagnetic superexchange of 120 meV along the c axis, the same size 
with the interactions in the cuprate superconductors, and ferromagnetic couplings smaller by an 
order of magnitude along a. CaIrC>3 can thus be regarded as a j = 1/2 quasi-one-dimensional 
antiferromagnet. While this qualitatively agrees with the stripy magnetic structure proposed by 
resonant x-ray diffraction, the detailed microscopic picture emerging from our study, in particular, 
the highly uneven admixture of ti g components, provides a clear prediction for resonant inelastic 
x-ray scattering experiments. 

PACS numbers: 71.15.Rf, 71.27.+a, 75.30.Et, 75.10.Dg 
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Iridium oxide compounds are at the heart of intensive 
experimental and theoretical investigations in solid state 
physics. The few different structural varieties displaying 
Ir ions in octahedral coordination and tetravalent 5d 5 va- 
lence states, in particular, have recently become a fertile 
ground for studies of new physics driven by the inter- 
play between strong spin-orbit interactions and electron 
correlations [1]. The spin-orbit couplings (SOC's) are 
sufficiently large to split the Ir t\ g manifold into well sep- 
arated j = 3/2 and j = 1/2 states. The former are fully 
occupied while the latter are only half filled. It appears 
that the widths of the j = 1/2 bands, w ~ Nt, where t 
is an effective Ir-Ir hopping integral and N is the coor- 
dination number, are comparable to the strength of the 
onsite Coulomb repulsion U. For smaller coordination 
numbers, such as N — 3 in the layered honeycomb sys- 
tems Li2lrC>3 and Na2lr03 and N = 4 for square-lattice 
Sr2lr04 and Ba2lr04 and the post-perovskite CalrO^, 
the undoped compounds are insulating at all tempera- 
tures and a j = 1/2 Mott-Hubbard type picture seems 
to be appropriate [1-3] . For larger coordination numbers 
like N = 6 in the pyrochlore structure, the t\ g Ir oxides 
are metallic at high T's and display transitions to antifer- 
romagnetic (AF) insulating states on cooling [4]. Here, 
a Slater type picture in which the gap opening is inti- 
mately related to the onset of AF order might be more 
appropriate, as proposed for example for Os oxide per- 
ovskite systems with JV = 6 [5]. Nevertheless, also more 
exotic ground states such as a Mott topological insulator 
[6] or a Weyl semimetal [7] have been proposed for the 
pyrochlore iridates. 

Besides the SOC's and Coulomb interactions, one im- 
portant additional energy scale is the magnitude of the 
splittings within the Ir t^g shell. Such splittings may 
arise from distortions of the O octahedron around the 



Ir site and also from the anisotropy of the farther sur- 
roundings. In this respect, highly anisotropic systems 
such as the post-perovskite Cali'03 provide an ideal play- 
ground for new insights into the interplay between SOC's 
and lattice distortions. In the layered post-perovskite 
structure, IrOg octahedra in the ac layers share corners 
along the crystallographic c axis, see Fig. 1, and have 
common edges along a [8]. The Ir-0 bonds along the 
corner-sharing chains are substantially shorter than for 
the edge-sharing links, 1.94 vs. 2.07 A, which suggests a 
stabilization of the xy-like orbital with respect to the xz 
and yz terms. Additionally, small splittings between the 
xz and yz components may in principle arise due to the 
a-b crystal anisotropy. 




FIG. 1. Local coordination and the main, 7r-type 5^2-2^3;- 
5d xz AF superexchange path for corner-sharing Ir06 octahe- 
dra along the c axis in CalrOa, see text. Nearest-neighbor 
O ions around the two Ir sites in the center of the figure are 
shown as small light-grey spheres. Ca ions above and below 
the Ir03 layer are not displayed. The O 2p tails of the active 
Ir d xz orbitals are overlapping at the bridging O site. 
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The precise nature of the relativistic 5c? 5 ground-state 
wave function obviously depends on the size of the Ir t2 g 
splittings. Further, the on-site 5d 5 electron configura- 
tion and the kind of Ir-O-Ir path along a given direction 
determine the sign and magnitude of the magnetic ex- 
change couplings. Ab initio calculations where all these 
effects and interactions are treated on equal footing are 
therefore highly desirable. For this, we focus on the post- 
perovskite phase of Calr03, which has received signif- 
icant attention from geologists in recent years because 
the same crystalline structure is adopted by MgSiOs in 
the lower Earth's mantle [9, 10]. In contrast to MgSiC>3, 
the post-perovskite phase of Calr03 can be synthesized 
even at ambient pressure. Since measurements are less 
complicated in ambient conditions, Calr0 3 has been used 
extensively as an analogue of MgSiC>3 post-perovskite for 
structural, mechanical, and transport studies [11-14]. To 
unveil the valence electronic structure of CalrOs we here 
carry out ab initio wave-function quantum chemical cal- 
culations. With no SOC's, splittings as large as 0.8 eV 
are found for the tig levels of CalrOa. The 5c? 5 rela- 
tivistic ground-state wave function therefore displays a 
highly uneven admixture of xy, xz, and yz character, 
which contradicts the interpretation of recent resonant 
x-ray diffraction (RXD) experiments [15]. For the corner- 
sharing octahedra, the quantum chemical calculations 
further predict large AF superexchange couplings [16] of 
120 meV, about the same size with the magnetic inter- 
actions in the parent compounds of the high-T c cuprate 
superconductors. The lowest d-d excited states are two 
doublets at approximately 0.7 and 1.4 cV, originating 
from the j = 3/2 quartet in an ideal cubic environment. 
This is the largest splitting reported so far for the j = 3/2 
quartet in d 5 systems and singles out CalrOs as an unique 
material in which both the SOC and the non-cubic crys- 
tal field splittings are of the order of half eV or larger. 

To investigate the local Ir c?-level splittings, multicon- 
figuration self-consistent-field (MCSCF) and multirefer- 
ence configuration- interaction (MRCI) calculations [17] 
are performed on embedded clusters made of one refer- 
ence Ir06 octahedron, ten adjacent Ca sites, and four 
nearest-neighbor (NN) Ir06 octahedra. That a realistic 
representation of ions at NN polyhcdra guarantees highly 
accurate results for the d-leve\ electronic structure at the 
central site was earlier pointed out for correlated 3d 1 , 3d 2 , 
and 3c? 9 oxide compounds [18-20]. The farther solid-state 
environment is modeled as a one-electron effective poten- 
tial which in an ionic picture reproduces the Madelung 
field in the cluster region. All calculations are performed 
with the MOLPRO quantum chemical software [21]. We 
employed energy-consistent relativistic pseudopotentials 
for Ir [22] and Ca [23] and Gaussian-type valence basis 
functions [22-25] (for details, see Supplementary Mate- 
rial [26]). For the ground-state calculations, the orbitals 
within each finite cluster are variationally optimized at 
the MCSCF level. All Ir ti g functions are included in the 



TABLE I. Ir t2 g splittings and relative energies of the spin- 
orbit tig states in Calr0 3 . MCSCF and MRCI results with 
and without SOC's for 5-plaquette clusters, see text. 



t'2 g splittings (eV) 



MCSCF 



MRCI 



-du 



txy dxz 



\j = l/2, mj = ±L/2) 
\j = 3/2, mj = i3/2) 
|j = 3/2,m 3 = ±l/2} 



0.63 
0.76 

0.00-0.12 
0.62-0.70 
1.25-1.27 



0.68 
0.83 

0.00-0.16 
0.66-0.76 
1.36-1.39 



active orbital space. Excitations within the Ir ti g shell 
are afterwards computed just for the central Ir06 octa- 
hedron while the occupation of the NN Ir valence shells 
is held frozen as in the MCSCF ground-state configura- 
tion. Spin-orbit interactions are further accounted for as 
described in Ref. [27]. The subsequent MRCI treatment 
includes all single and double excitations from the O 2p 
orbitals at the central octahedron and the Ir 5c? func- 
tions. To partition the O 2p valence orbitals into two 
different groups, i.e., at sites of the central octahedron 
and at NN octahedra, we employ the orbital localization 
module available with molpro. Crystallographic data 
as reported by Rodi and Babel [8] is used. 

With no SOC's, we find Ir d xy —d yz and d xy —d xz split- 
tings of 0.63 and 0.76 eV by MCSCF calculations and 
slightly larger values of 0.68 and 0.83 eV at the MRCI 
level, with the xy level the lowest (see the first lines in Ta- 
ble I). The j = l/2 ground-state wave function therefore 
has dominant xz character, as much as 75% by MRCI 
calculations with SOC (MRCI+SOC). The four compo- 
nents of the j = 3/2 quartet are split into groups of two 
by 0.6-0.7 eV (lower lines of Table I). The lower two 
j = 3/2 terms have 72% yz character, the highest j = 3/2 
components have 90% xy character. 

When SOC's are not accounted for, the coupling to 
spin moments at NN Ir sites gives rise for each partic- 
ular orbital occupation to one sextet, four quartet, and 
five doublet states. All those configuration state func- 
tions enter the spin-orbit calculations. 32 different spin- 
orbit states are therefore generated for each j configura- 
tion at the central Ir site and this causes a finite energy 
spread in Table I. A simpler and more transparent pic- 
ture can be obtained by replacing the four Ir 4+ d 5 NN's 
by closed-shell Pt 4+ d 6 ions. In that case, the MRCI ex- 
citation energies for the j = 3/2 states at the central Ir 
site are 0.67 and 1.38 eV. An instructive computational 
experiment which we can further make is to replace the 
distorted post-perovskite structure by a hypothetical cu- 
bic perovskite lattice in which the Ir-0 bond length is 
the average of the three different bond lengths in post- 
perovskite CalrOa, i.e., 2.03 A. This allows to directly 
extract the SOC A for the Ir 4+ 5c? 5 ion because the split- 
ting between the spin-orbit j = l/2 doublet and the four- 
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TABLE II. NN effective magnetic couplings in CaIrC>3, see 
text. Positive values denote AF exchange. 



J (meV) 


MCSCF 


MRCI 


MRCI+SOC 


Jc 


91.2 


163.4 


121.0 


Ja 


-2.3 


4.2 


-7.3 



fold degenerate j = 3/2 state is 3A/2 in cubic symmetry 
[28]. The doublet-quartet splitting comes out as 0.70 eV 
by MRCI+SOC calculations for the idealized cubic per- 
ovskite structure. This yields a spin-orbit coupling con- 
stant A = 0.47, in agreement with values of 0.39-0.49 eV 
earlier extracted for Ir 4+ impurities from electron spin 
resonance measurements [29]. 

The results of our ab initio calculations, i.e., the un- 
usually large t 2g splittings and the highly uneven admix- 
ture of xz, yz, and xy character for the spin-orbit t\ g 
states, are in sharp contrast to the interpretation of re- 
cent RXD experiments on CalrOa [15]. In particular, 
the authors of Ref. [15] conclude from the RXD data 
that the deviations from an isotropic xz-yz-xy picture 
are marginal, as also found for instance in the layered 
square-lattice material Sr2ir04 [1]. While for Sr2ir04 
we also obtain [30] rather small t 2g splittings of about 
0.1 eV, comparable weights of the xz, yz, xy compo- 
nents in the j — 1/2 ground-state wave function, and NN 
superexchange interactions in good agreement with reso- 
nant inelastic x-ray scattering (RIXS) data, for the post- 
perovskite CalrOs the quantum chemical results point 
in a different direction. These conflicting findings ur- 
gently call for Ir L-edge RIXS measurements on Calr03. 
The RIXS experiments have recently emerged as a most 
reliable technique for exploring the charge, spin, and or- 
bital degrees of freedom of correlated electrons in solids 
[31, 32]. 

The presence of both corner-sharing and edge-sharing 
octahedra and the large splittings within the t 2g shell 
have important implications on the nature and the mag- 
nitude of the magnetic interactions. To determine the 
latter, we designed 8-octahedra clusters including two 
active Ir 4+ d 5 sites. The six 5d NN's were modeled as 
closed-shell Pt 4+ d e ions. All possible occupations were 




FIG. 2. Sketch of the NN magnetic interactions in CaIrC>3. 



allowed within the set of t 2g orbitals at the two active Ir 
sites in the MCSCF calculations, which gives rise to nine 
singlet and nine triplet states with no SOC. The MCSCF 
wave functions were expressed in terms of orbitals opti- 
mized for an average of those singlet and triplet states. 
All eighteen singlet and triplet states entered the spin- 
orbit calculations, both at the MCSCF and MRCI levels. 
In MRCI, single and double excitations from the Ir ti g 
orbitals and the 2p orbitals of the bridging ligand were 
included. Such a computational scheme provides J val- 
ues in very good agreement with RIXS measurements in 
the layered iridate Sr 2 Ir04 [30]. 

Remarkably, with no SOC's, the NN AF exchange con- 
stant J c for octahedra sharing corners along the c axis 
is as large as 91.2 meV by MCSCF calculations and 
163.4 meV by MRCI, see Table II. The t 2g hole or- 
bitals are here the xz components, as discussed above. 
Although the IrOe octahedra are tilted somewhat due to 
rotations about the a axis, see Fig. 1, this does not af- 
fect much the 7r-type 5d xz -2p x -5d xz overlap and superex- 
change [16] interactions. An AF J of 163.4 meV is even 
larger than the values for cr-type superexchange paths in 
two-dimensional S — 1/2 3d Cu oxide superconductors 
[33] and points to the role of the larger spatial extent of 
5d functions in iridates. With SOC's accounted for, the 
effective AF coupling J c is reduced to 121.0 meV (last 
column in Table II). The spin-orbit ground-state wave 
function now acquires sizable yz character and since the 
tilting of the IrOg octahedra does affect the out-of-plane 
5d yz -2py-5d yz superexchange path, the effective magnetic 
coupling between the j = 1 / 2 sites is smaller. 

For edge-sharing octahedra along the a axis, we find a 
weak NN ferromagnetic (FM) coupling J a sa —7.3 meV in 
the spin-orbit MRCI calculations. The three components 
of the jtot = 1 triplet state in the two-site problem display 
small splittings of « 1 meV (see Supplementary Mate- 
rial [26]), which shows that anisotropic non-Heisenbcrg 
terms [34] should also be considered when mapping the 
ab initio data onto an effective spin Hamiltonian. In esti- 
mating the Heisenberg coupling constant J a we neglected 
such terms and used an average of the energies of the 
jtot = l triplet components, JmE(j tot = 1) -E(j tot =0). 
The ab initio results, i.e., strong AF interactions along c 
and weak FM couplings along a (see Fig. 2), are in full 
agreement with the observed striped AF magnetic order- 
ing and very large Curie- Weiss temperature [15]. The 
MRCI+SOC J values thus provide a firm quantitative 
basis in understanding the magnetic properties of this 
material. 

We have employed, in summary, a set of many-body 
quantum chemical calculations to unravel the 5c? elec- 
tronic structure of the post-perovskite iridate CalrOs. 
While electronic structure calculations have been earlier 
performed within the local density approximation to den- 
sity functional theory and the role of electron correla- 
tion effects has been pointed out [35] , we here make clear 
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quantitative predictions for the Ir t2 g splittings, character 
of the entangled spin-orbit wave function, and magnitude 
of the NN superexchange interactions. Our results single 
out CalrOa as a 5c? 5 system in which lattice distortions 
and spin-orbit couplings compete on the same energy 
scale to give rise to highly anisotropic electronic structure 
and magnetic correlations. In particular, the present ab 
initio investigation yields a different picture than previ- 
ously derived on the basis of resonant x-ray diffraction 
experiments for the relativistic ground-state wave func- 
tion. The large t 2g splittings that we have found from 
the calculations give rise to dominant xz hole character 
and remarkably strong AF interactions along the c axis. 
Further, the large AF couplings along c and weak FM 
exchange along a characterize Calr03 as a j = l/2 quasi- 
one-dimensional antiferromagnet. 
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SUPPLEMENTARY MATERIAL 
COMPUTATIONAL DETAILS 
Ir d-level splittings 

To determine the splittings within the Ir i 2 g levels and 
spin-orbit t\ g states we employed large 5-octahedra clus- 
ters embedded in arrays of point charges fitted to re- 
produce the crystal Madelung field in the cluster region. 
Quadruple-zeta basis sets from the molpro library were 
applied for the valence shells of the central Ir ion [1] and 
triple-zeta basis sets for the ligands [2] of the central oc- 
tahedron and the nearest-neigbor Ir sites [1]. For the 
central Ir ion we also used two polarization / functions 
[1]. For farther ligands in our clusters we applied min- 
imal atomic-natural-orbital basis sets [3]. All occupied 
shells of the Ca 2+ ions were incorporated in the large- 
core pseudopotentials and each of the Ca 4s orbitals was 
described by a single contracted Gaussian function [4] . 



TABLE III. Relative energies (meV) of singlet and triplet 
states for two adjacent Ir sites for corner-sharing octahedra. 
The singlet state is taken in each case as reference. 





MCSCF 


MRCI 


MRCI+SOC 


Singlet 


0.0 


0.0 


0.0 


Triplet 


91.2 


163.4 


120.0 








120.2 








122.8 



TABLE IV. Relative energies (meV) of singlet and triplet 
states for two adjacent Ir sites for edge-sharing octahedra. 
The triplet or lowest triplet component is taken in each case 
as reference. 





MCSCF 


MRCI 


MRCI+SOC 


Triplet 


0.0 


0.0 


0.0 








0.4 








1.5 


Singlet 


2.3 


-4.2 


7.9 



Effective magnetic couplings 

Calculations for the isotropic magnetic exchange inter- 
actions were performed on 8-octahedra clusters with two 
active Ir 4+ 5cZ 5 sites. To make the analysis of the inter- 
site spin couplings straightforward, the six Ir 4+ t\ g NN's 
were replaced by closed-shell Pt 4+ t% g species. The basis 
functions used here were as described above. The only 
exception was the O ligand bridging the two magneti- 
cally active Ir sites, for which we employed quintuplc- 
zeta valence basis sets and four polarization d func- 
tions [2] . Relative energies for the lowest 5 to t = and 
jtot = singlet and S t ot = 1 and j to t = 1 triplet states for 
corner-sharing and edge-sharing octahedra are listed in 
Table III and Table IV, respectively. In the spin-orbit 
calculations the triplet components display small split- 
tings that indicate the presence of small non-Hciscnbcrg 
terms. In deriving the effective Heisenberg coupling con- 
stants J c and J a we neglected such terms and used an 
average of the energies of the jtot = 1 triplet components, 
J~E(j tot = l)-E(j tot = 0). 
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